Stochastic origin of Gompertzian growths 
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This work faces the problem of the origin of the logarithmic character of the Gompertzian growth. We 
show that the macroscopic, deterministic Gompertz equation describes the evolution from the initial state to the 
final stationary value of the median of a log-normally distributed, stochastic process. Moreover, by exploiting a 
stochastic variational principle, we account for self-regulating feature of Gompertzian growths provided by self- 
consistent feedback of relative density variations. This well defined conceptual framework shows its usefulness 
by allowing a reliable control of the growth by external actions. 

PACS numbers: 



I. INTRODUCTION 

The Gompertz model, in its original conception, was born 
as phenomenological one namely describing the observed age 
tables of humans |[l2tl . In fact, B. Gompertz concluded his em- 
pirical studies of tables introducing the distribution of human 
ages for a given community, the now well known function 



P(r) 



(1) 



where a > 0, (3 < and c are constant. 

It is interesting to note that Gompertz posed at the core of 
his deduction the properties of geometrical progression; "This 
law of geometrical progression pervades, in an approximate 
degree, large portions of different tables of mortality". The 
relevance of the geometrical progression in the framework of 
the natural phenomena in a variety of experimental environ- 
ments was pointed out at the end of nineteenth century in the 
works of Galton U3] and McAlister They showed that 
the geometrical mean (median) describes the behavior of a 
large set of natural phenomena better than the arithmetic one. 
We note that Gompertz law according to the Gompertz obser- 
vation emerges from the equilibrium between geometrical se- 
ries associated to degradation and an arithmetical progression 
ruling indefinite (Malthusian) growth having the experimental 
observations to be made on suitable time intervals. The char- 
acteristics of simplicity of Gompertz law due to this general 
and profound mathematical framework, attract the attention 
of nascent biological disciplines, where the growth studies go 
back to the thirty years of last century. Since this distribution 
has had so remarkable success in a variety of very different 
situations that a lot of literature refers it simply as "law of 
growth". 

The laws of growth of natural systems, and the deep ori- 
gin of their characteristic scales of, e. g., length, mass, en- 
ergy, or numerosity, are, now, intensively investigated in many 
branches of science, such as biomedicine (a 1 1 €41 - economy 
1 9], population dynamics [18], astrophysics and cosmology 
\M- So, starting from 1930 the Gompertz equation has be- 
came one of the most used tool to account for mechanisms of 
growth in a variety of systems in many fields fl-H HI [l8tl . If 
we exclude, for example, from incomes distribution the little 



percentage of rich people ( 1 %) Gompertz law is the fitting dis- 
tribution. Obviously so general applicability has developed a 
very interesting debate regarding the origin of its logarithmic 
structure. The arguments called down to this aim look at the 
main aspects of the underlying systems i.e., biological social 
or/and economics. 

General theory of dynamical systems, quiescence, cell ki- 
netics theory, entropic and thermodynamical arguments have 
been advocated and illustrated by many authors in a variety of 
interesting papers along many years introducing also suitable 
generalizations and connections with other growth models as 
logistic one. A good synthetic description with a large bibli- 
ography can be considered that of Bajzer et al (1997)]2|]. 

We remark three main aspects of the delineated problem 
that we consider relevant. 

The first one; although, it can be significant to start attach- 
ing the problem from specificity of a discipline, in fact this can 
enlighten nodal points, the arguments bringing to logarithmic 
behavior must be so general as well as is the applicability of 
Gompertz law; 

The second one; Gompertz curve cannot be other than a 
"suitable" description of a mean behavior of systems under 
studies that are all characterized by a basic stochasticity. 

The third one; many of the considered systems reach the 
limiting size exploiting self-controlled evolutions. We take 
these considerations as starting points to recovering the fea- 
tures, in particular logarithmic behavior, of Gompertz equa- 
tion. Before to go into the details of our deduction it is useful- 
ness to make some preliminary consideration and as first step 
to give a description of deterministic Gompertz equation also 
to establish our notation. 



II. GOMPERTZ EQUATION 

The standard form of the deterministic Gompertz equation 

is: 



dz 
~dl 



z 



(2) 



where z describes the "size" of some quantity characteriz- 
ing the system, f3 and a denote two positive constants with the 



2 



dimensions of the inverse of time, and z is a constant which 
the same dimensions of z. 
Eq. (ffj) can be recast as: 

^—L = -a Ins, (3) 
at 

where s(t) = z(t)/z ao , and z x = zexp(j3/a). The Gom- 
pertz equation is then associated to four parameters (all depen- 
dent on the specific system): a, (3, z, and the initial condition 
("scale") z(0) = Zq. Its solution is: 

z(t) = z 00 cxp[(ln 7 )-e- Q *], (4) 

where 7 = Zq/ 1 z^. It is immediately verified that this solu- 
tion always approaches monotonically in time z^: depending 
on the conditions zq < z^ = zcxp/3/a (7 < 1), or 
zq > Zoo = zexp(3/a (7 > 1), the system monotonically 
grows or monotonically decreases, respectively, from the di- 
mension z to the dimension z^, approaching the asymptotic 
value Zoo with the characteristic time a -1 . It is worth notic- 
ing that the solutions of the eq. (01 satisfy some, somewhat 
simple and self evident properties, deriving by peculiar fea- 
tures of logarithmic function. These properties justify why the 
Gompertz equation plays a major role among all the equations 
describing growth phenomena. 

If we characterize a solution s(t) by the pair of values of its 
parameters (7, a), the properties are the following: 

1. the product si(t) ■ S2(t) of two solutions with parame- 
ters (71, a) and (72, a), respectively, is again a solution 
with parameters (71 • 72, a), 

2. if s(t) is a solution with parameters (7, a), then s a (t), 
with a € R is a solution as well, with parameters 

(7 Q ,a). 

3. the constant function s(t) = 1 is a (trivial) solution. 

Note that, for a = —1 we obtain that the inverse of s(t), 
s _1 (f) is a solution associated to (7 -1 , a). Then it is straight- 
forward to verify that if for example the original solution 
is obtained with zo > z^o (aggregation process) s _1 (i) de- 
scribes the time-reversed (fragmentation) process from z^ to 
Zq. Note also that this property, together with properties 1), 
3) makes the set of solutions of the Gompertz equation an 
Abelian group. Note finally that the property 3) implies that 
any other quantity M(t) linked to z(t) by an allometric rela- 
tion (M (t) ~b ■ z a (t)) satisfies the same Gompertz equation 
with modified parameter 7° 



view was that diseases and disfunctions accumulate slowly 
along the time damaging multiplicatively the human bodies. 
The extensive analysis of Jones showed that Gompertz equa- 
tion applies exactly to people that have not eliminated the 
first cause of diseases, i.e. hygienic condition. A good mean 
improvement of these last one, namely slight, modifies the 
general behavior. Then, following Jones, at the basis of the 
analysis of life expectancy, there is a stochastic process built 
with independent random variables (diseases and/or social and 
economic condition) that add multiplicatively. Moreover at 
the end of sixties years of the last century a detailed statis- 
tical analysis performed by Sachs showed that physiological 
parameters like blood pressure, tolerability of medicaments, 
body size survival rate are lognormally distributed l2lll . Fi- 
nally it is worthwhile to note that the lognormal i.e. geo- 
metric Brownian motion appears to be asymptotic to a vari- 
ety of branching processes introduced to describe cell systems 
growth 011. 

On the other hand it is well known that lognormal dis- 
tribution arises in a variety of classification procedures and 
in physical and biological systems when natural genesis in- 
volves repeated breakages or aggregations. Very relevant, in 
1941, Kolmogorov lfl5f1 has shown that when the frequency 
of aggregation-disaggregation in a growth process is indepen- 
dent of the size of the constituents, the asymptotic size dis- 
tribution of the aggregate should tend to be lognormal. This 
so general multiplicative behavior of basic randomness of un- 
derlying processes has been indicated as a "Multiplicative 
Gestaltungs-Principle of Nature". It is also assumed as its 
possible accounting a general property of "coherence" of nat- 
ural systems l22tl . The simple considerations outlined above, 
at the light of the two first points of our introductory section 
namely basic stochasticity of systems and large applicability 
of the Gompertz model, bring us to postulate that the Gom- 
pertz equation is the deterministic one emerging from geo- 
metric Brownian motion, that is the general stochastic struc- 
ture associated with a variety of systems. The necessary step 
at this point is to give a simple mathematical procedure con- 
necting stochastic structure with deterministic one. 

B. Mathematical procedure 

The main feature of lognormal can be so summarized: a 
random variable X is said to have a lognormal distribution 
with suitable parameters associated to mean and variance if 
ln{(X (t))} is normally distributed. Consider then the diffu- 
sion process X(t) taking non negative values, and satisfying 
the Ito differential equation: 



III. GOMPERTZ EQUATION AS EVOLUTION EQUATION 
FOR THE MEDIAN OF GEOMETRICAL BROWNIAN 
MOTION. 

A. Lognormality as basic in a variety of natural systems 

In 1947 H.R. Jones linked the problem of mortality to life 
expectancy and ageing processes fTill . His germinal point of 



dX(t) = {[-X(t)-aX(t) In [-Ll]}dt+y/i;X(t)dW(t), 

(5) 

where a and v are positive constants, K has the same dimen- 
sions of X(t), and dW(t) denotes a Gaussian stochastic pro- 
cess (Wiener process) with zero mean and variance dt. Here 
X(t) can describe any quantity. X(t) is a multiplicativy diffu- 
sive process (geometrical brownian motion) and its Ito equa- 
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tion can be recast as: 



-X(t). 



dX(t) = -aX(t) In dt + s/vX(t)dW(t), (6) 



where 



K = Ke~° 



We now prove that this process generates the deterministic 
Gompertz equation. Let us in fact associate to X(t) a new 
process Y (t) defined by the relation 



X(t) 



(7) 



where Y(t) is adimentional and takes values ranging on 
( — oo, +oo) when X(t) takes values on [0, oo ). By exploit- 
ing Eq. (O and definition|7] we can compute c?(exp Y(i)): by 
Ito's lemma obtaining: 



dY(t) = -aY(t)dt + y/vdW(t), 



(8) 



We see that Y(t) is an Omstein-Uhlenbeck process . Its prob- 
ability density p(y,t;y ,Q) = p(y,i), satisfies the Fokker- 
Plank equation 

dtp(y, t) = ad y [yp{y, t)] + ^v r d 2 y p(y, t), (9) 
and can be exactly computed for any initial condition 1 1 1]: 



due- my xo{ue- at )e 



(10) 



where Xo( u ) i s tne characteristic function of the initial proba- 
bility. The moments < Y n (t) > of the process Y(t): 



we see that z(t), which has the same dimensions of X{t), sat- 
isfies the Gompertz equation, where the constant K is iden- 
tified with Zqo, i- e. with the asymptotic value of z(t). But 
what about the meaning of z(f) ?. We nown provide a precise 
meaning to z(t). Let us suppose that the initial probability of 
the process Y(t) is Gaussian (for example, the solution is the 
fundamental one with an initial delta function condition). By 
Eq. ( [Tot , the solution at any time is Gaussian: 

p(y, t)dy = J— e-^^^dy, (16) 

where a(t) =< [Y(t) - fi(t)] 2 >, and p,(t) =< Y(t) >. By 
exploiting the relation ©, we can recast Eq. ( fTTI ). obtaining 
the probability of the process X(t) as 

p( x ,t)dx = - 1 ,[-^)0"(^)^w) 2 ] (17) 

where, by solving the first two branching equation (fTSl i). mean 
fj,(t) and variance cr(t) are given by: 



H(t) = pL e 



<?{t) = <J e + — {1- e ) 
la 



(18) 



(19) 



<Y n (t)> = / dyy n p{y,t), (11) 

J — oo 

satisfy the set of branching equations 



— < Y n (t) >= -na < Y n {t) > +^a r n(n-l)a < Y n ~ 2 {t) 

(12) 

which can be solved, iteratively at any finite order. We now 
focus on the case n = 1: 



dt 



< Y(t) >= -a < Y(t) >, 



(13) 



exploit the relation (0, and obtain 



If we define the time-dependent quantity z(t) by the relation 

1 <'4^ i ]>, (15) 



K 



exp 



We see that the process X{t) /K is lognormally distributed. 
But it is well known that in a lognormally distributed process 
the mean of the logarithm of the process is the logarithm of 
the median of the process. Then, we can conclude that if the 
process Y{t) is initially Gaussian, i.e. if the process x(t)/K 
is initially lognormally distributed, this last process remains 
lognormally distributed at any time, and the variable z(t)/K 
(which we denoted z(t) / z^) is the median of this process. 

In conclusion we have proved that: the deterministic Gom- 
pertz equation is the macroscopic consequence of a lognor- 
mally distributed, diffusion process X(t); the macroscopical 
size z( t) whose evolution is ruled by the Gompertz equation is 
the median of process X(t). 

It follows that the multiplicative stochastic process (a stan- 
dard geometric brownian motion) that many times is intro- 
duced as "Gompertz stochastic tumor growth model" is not 
>due to an extra noise disturbing the Gompertz growth, but it is 
itself the origin of the deterministic growth. 

To further support this conclusion we remark that deter- 
ministic Gompertz model was extracted by B. Gompertz, 
just looking to the properties of a geometric progression that 
emerges from the mortality tables; i.e. the Gompertz func- 
tion has its natural interpretation as median of a multiplicative 
process. 



C. Short remark on the observability of microscopic 
parameters 

The underlying process X(t) is characterized by three pa- 
rameters: a, v, K. Their determination completely define the 
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process. As a first observation, we note that the Malthusian 
parameter j3 in the standard form (O of the Gompertz equa- 
tion is provided, in a consistent way, by the diffusion param- 
eter: f3 = v/2. We remark also that the drift parameter a 
is clearly macroscopically observable by fitting the macro- 
scopic size growth of the selected systems. Therefore, if one 
is able to provide a method to extract by observational data 
the diffusion parameter v, the whole process X(t) can be re- 
constructed. To this end, we consider that the observed size 
(i. e., the median of the process) must displace a range of 
variability; in fact 

\nz{t)/ ZoQ = \nz{t)/K =< \nX(t)/K >=< Y(t) > 

(20) 

is the mean of a Gaussian process with variance (width) cr(t). 
And, in particular, when the last stage K = is reached, the 
width (i.e variability of the size) becomes that of the stationary 
distribution (see Eq. (fT9]l). If one performs a statistics of the 
observed sizes, one could find that the sizes of a system at the 
last stage range from a minimum one z m %n to a maximum one 
Zmax, with the following relations: 



^oo V^min ' Zmax 7 (21) 

(just geometric mean), 

r 

and 

, Z max /TON 

^ = aln . (23) 

Zjnin 

Summing up: a, z max , z m i n can be extracted by observational 
data, and, in principle, their values provide an estimate of all 
the other parameters of the underlying process; in particular, 
of the order of magnitude of the diffusion parameter v. 



IV. CHARACTERIZING THE GOMPERTZIAN GROWTH 
PROCESS 

The conclusions reached in the previous section ascribe the 
Gompertz equation within the framework of "Multiplicative 
Gestaltung-Principle of Nature", i.e. to the ubiquitous na- 
ture of lognormal distribution, which characterizes a variety 
of natural system, i.e. those in which emerge characteris- 
tic scales of "coherence". The "actual" Gompertz equation 
is then the stochastic one. The stochastic description how- 
ever does not exhaust all the aspects of Gompertzian model, 
namely stochasticity allows spontaneous growth until a final 
size if and when systems undergo to a stationary state. The 
Gompertzian growth is characterized, on the contrary, as re- 
marked in the third point of introductory section, by a self- 
controlled evolution ruled by variation of density. Now we 
incorporate in our description this relevant behavior of grow- 
ing systems. 



A. Dynamical updating of geometrical Brownian motion 

In standard treatment the drift terms in a Fokker-Planck 
equation is a function given a priori; backward and forward 
evolutions are described with two different equations. The 
systems spontaneously approach a stationary state when bal- 
ancing is reached between the stochastic term (Wiener pro- 
cess) and deterministic (drift); consequently there is not, for 
example, time reversal invariance. In order to have self- 
control (feedback) we must to add to Fokker-Planck equation 
a dynamical equation updating the drift and allowing time re- 
versal behavior. The equation is the stochastic equivalent of 
F = ma, and it can be written also as the equation describing 
a suitable interface l20ll . 

We remark, that, being our systems self-controlled the dy- 
namical update must contain self coupling non linearity. We 
use a stochastic variational principle theory. We recall this 
theory briefly: deterministic dynamic evolutions are char- 
acterized by two independent principles, the first kinematic 
the second dynamic. The kinematic principle is provided by 
standard differential rules, and the dynamic one by a varia- 
tional principle, i.e. the Lagrangian principle. The theory 
of stochastic variational principle assumes Ito's equation as a 
kinematic rule and the Lagrangian variational principle as dy- 
namic rule. The variation is made by considering conditional 
expectation. As a consequence: the configuration of our sys- 
tems is described, in general, by a vectorial Markov process 
taking values in 5ft 3 . This process is characterized by a 
probability density p(r, t) and a transition probability density 
p(r, t\ r', i'), and its components satisfy an Ito stochastic dif- 
ferential equation of the form 

d^(t) = v {+)j ^(t),t)dt + dr)j{t), (24) 

where are the components of the forward velocity field. 
As already observed here the fields i>(+)j must not be given a 
priori, but play the role of dynamical variables and are con- 
sequently determined by imposing a specific dynamics. The 
noise rj (t) is a standard Wiener process, D is the diffusion 
coefficient. We indicate by E t the conditional expectations 
with respect to In what follows, for sake of notational 
simplicity, we will limit ourselves to the case of one dimen- 
sional trajectories, but the results that will be obtained can be 
immediately generalized to any number of dimensions. We 
will suppose for the time being that the forces will be defined 
by means of purely configurational, possibly time-dependent 
V (x, t), potentials, this includes also a non linear potential 
functional of the density of the process. A suitable definition 
of the Lagrangian and of the stochastic action functional for 
the system described by the dynamical variables p and 
allows to select, the processes which reproduce the correct 
dynamics JfJ, [TH [l^] . In fact, while the probability density 
p(x, t) satisfies, as usual, the forward Fokker-Planck equation 
associated to the stochastic differential equation 

d t p = Dd 2 x p - d x {v {+) p) = d x {Dd x p - v (+) p) (25) 

the following choice for the Lagrangian field 
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in 

L(x, t) = —vf +) (x, t) + mDd x v (+) (x, t) - V(x, t) (26) 



enables to define a stochastic action functional 

rti 



A 



Et[L(£(t),t)]dt 



(27) 



which leads, through the stationarity condition 5 A = 0, to the 
equation 



= 0. 



(28) 



2m yfp 
The field S(x, t) is defined as 

S(x, t) = - J' E [L(S(s),8) | m = x] ds + 

+E[S 1 (£(t 1 ))\£(t)=x] (29) 

where Si(-) = S( ■ ,ti) is an arbitrary final condition. 
This equation is the well know equation of interfaces theory 
where the term depending on density p represents the con- 
tribute due to surface tension. 

By introducing the function R(x, t) = p(x, t) and the de 
Broglie ansatz 



i>(x,t) = R(x,t)e zS{x ^ /2mD 
equation d28l i takes the form 



dt S + ^ + V-2mD* 9 -^=0 



2m 



R 



(30) 



(31) 



and the pair of real equations d25l l and ( |3"TT l are equivalent to 
the single linear equation for tjj 



i{2mD)d t ip = - 2mD 2 d 2 ?p + Vip, 



(32) 



This connects our dynamic equations with an ordinary 
eigenvalue problem on Hilbert space. 

Note that the observables are the density p, which here rep- 
resents just the (mass) density of the system, and the drift 
mean velocity of the system v; where the connection with 
the pair i)),v is provided, at every point and at every time, 
by p = \ij)\ 2 , and v = d x S/m. Note also that 



-'(+) 



d x S + mD 



d x R 
R ' 



(33) 



If we choose the potential in the form of V(x) = f(p), where 
f(p)) is a nonlinear functional of density (e.g. ln(p)), we 
obtain a dynamic system that provides self-regulation by the 
self consistent feedback action of relative density variations 
typically associated to osmotic phenomena across interfaces. 
So well this very general model contains the third of relevant 
aspects considered to describe growing systems, in particular 
cells proliferation in solid tumor. Lognormal stochastic back- 
ground, self-control and interface theory are the constitutive 



elements of this model, which leads to the deterministic Gom- 
pertz equation as equation for the median. This approach, be- 
ing conceptually well-founded, allows to face, in principle, 
the growth problems as a controlled one namely it is possible 
to intervene with outdoor control. For example our method 
poses the reduction of tumoral mass from a new and interest- 
ing point of view as we describe in the next section. 



B. Controlled growth 

In this section we move on to implement the controlled evo- 
lution. In fact we exploit the transition probabilities of the 
Gompertz self-controlled evolution to model controlled evo- 
lutions from a given initial state to arbitrarily assigned final 
states. We start by observing that to every solution ip(x,t) 
of the Fokker-Planck equation ( 1251 . with a given V(+)(x,t) 
and constant diffusion coefficient D, we can always associate 
a "wavefunction" descriptions of eq.(l32l. of a dynamical sys- 
tem. To this aim, it is sufficient to introduce a suitable time- 
dependent potential V c (x,t), by exploiting the wave equation 
(l32l as a control equation Q. 

Here we quickly recall those elements of the controlling 
procedure which are needed for our aim, referring to Cufaro 
et al. 1999 @] for further details. Let us consider a solution 
p(x, t) of the Fokker-Planck equation, with a given (x, t) 
and a constant diffusion coefficient D; let us introduce the 
functions R(x, t) and W(x, t) defined by 



p{x,t) = R 2 (x,t) , 



and remind that the relation 



v (+) (x,t)=d x W(x,t), (34) 



d x S + mD 



d x R 
R 



d x p 



= d x S + mD/2—=d x (S + mD/2\np) (35) 
P 

must hold, where p is an adimensional function (argument of 
a logarithm) obtained from the probability density p by means 
of a suitable and arbitrary dimensional multiplicative constant. 
If we now impose that the function S(x, t) must be the phase 
of a wave function, we immediately obtain from the Eqs. (l34b 
and I 



S(x,t) = mW{x,t) - ^ In -6(t), (36) 

which allows to determine S from p and vr + \ (namely W) up 
to an additive arbitrary function of time 0(t). However, we 
must ensure that the wave function (f34t with R and S given 
above, is a solution of the wave equation (1321 . Since S and R 
are now fixed, equation (l32l must be considered as a relation 
(constraint) defining the controlling potential V c , which, after 
straightforward calculations, yields 



+V( +) d x In p) 



V c (x,t) = mD 2 d 2 x \np + mD{d t \np 
mv 2 



(+) 



md t W + 9. 



(37) 
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Of course, if we start with a wave function ip(x,t) associ- 
ated to a given time-independent potential V(x), the self- 
consistency is ensured, and this formula always yields back 
the given potential, as it should. 

We now leave the general way, and focus on the very in- 
teresting case, useful for our goal, of simple controlling po- 
tentials able to produce a evolution which can vary, and in 
particular reverse, the growth trend. We start by observing 
that, to our purposes, it is expedient to handle the Gaussian 
Ornstein-Uhlenbeck process Y{t), Eq. ([§), because this is a 
simple task, and it is equivalent to modify the underlying non 
Gaussian process X(t). In fact, being < Y(t) >= fi(t) =< 
In (X(t)/K) >= lnz(t), where z(t) is the median describ- 
ing the macroscopic size, the monotonicity of the logarithmic 
function ensures that a reduction of fj,(t) leads to a correspond- 
ing (multiplicative) reduction of the size. The probability dis- 
tribution for Y(t) is the Gaussian ( TP71 ). characterized by the 
two time-dependent parameters <r(t)). Let us suppose 

now that the time evolution at some instant t has led to the 
pair of values /j.(t) = (j,q, a(t) = <jq, and that we aim to re- 
duce this values in such a way that, after some characteristic 
time t, they becomes \i\ < [Iq, a\ < <7o. As proved in 
Cufaro et al.Q], in this Gaussian instance the controlling po- 
tential to be applied is harmonic, and has the form: 

V c (x, t) = Y [w 2 {t)x 2 - 2a(t)x + c(t)] . (38) 

This harmonic potential is completely determined by the 
choice of mean and variance, fj, c (t),a c (t) , of the Gaussian 
process self consistently generated by the controlling har- 
monic potential d38l . In fact, the time-dependent coefficients 
ui 2 (t) , a(t) , c(t) are all functions of /i c (t),a c (t) and of their 
time derivatives, but a further free function of time that can be 
exploited to simplify the expression of the potential or of the 
phase. Being we here merely interested to outline the general 
conceptual frame allowing control, we avoid burdening this 
subsection with the explicit expressions of the time-dependent 
coefficients, which can be found in Cufaro et al. (g). We only 
sum up again the procedure. One chooses the form of the con- 
trolling mean and variance /J, c (t), <r c (t), which in the charac- 
teristic time r go to the final reduced values n\,(T\\ then one 
inserts their expressions, and a suitably chosen form of the 
further free function of time, in the time dependent param- 
eters uj 2 (t), a(t), c(t), so obtaining an harmonic controlling 
potential that, applied to the system, drives it towards the re- 
duced mean and variance. We comment here on an important 
aspect: one can choose many functions fi c (t),a c (t) leading 
to the reduced mean and variance. If we set fi c (t) = fiof(t), 
and a c (t) = aog(t), the functions f(t),g(t) must satisfy the 
constraints: 

/(°) = !> /(°°) = Mo/^i 5 .9(0) = 1, g(oo) = cro/a!. 
A possible form for g(t), for example, could be: 

1 + 6exp (i/r) 
c + dcxp (t/r) ' 



with c + d = 1 + b: d = b- (erg / 'fi ) ■ Another possible choice 
can be found in [g] . The important point is that the choice of 
this functions, together with that of the further free function 
of time, must be made in such a way to define a controlling 
potential V c (x, t) which can be effectively engineered and ap- 
plied to the system. It is also evident that, from the point of 
view of practical implementations, one can resort to suitable 
approximations which can anyhow realize the goal within a 
permissible error. 



V. CONCLUSION 

The ubiquity of the Gompertz equation, raises a very in- 
teresting question in the framework of growth phenomena, 
namely the key presence of the logarithmic function ruling 
the nonlinear growth. In this work, we have shown that the 
ubiquity of this equation and its logarithmic regulation are the 
macroscopic expression of the ubiquity of the log-normal dis- 
tribution. In fact, we have proved that Gompertzian growths 
are generated by a log-normally distributed stochastic pro- 
cess, being the macroscopic Gompertz equation the evolu- 
tion in time of the median of the process. The median then 
describes the macroscopic size of the growing system. We 
remark, therefore, that the growth is not, as often supposed, 
disturbed by a stochastic noise, but that the stochastic pro- 
cess is its origin, and its guidance at any time. This scheme 
agrees with the claim of Galton ifToh and McAlister Q3| who 
at the end of the nineteenth century have shown that many 
natural systems are well described by log-normality; their ac- 
tual behavior is then described by the median rather than the 
mean, thus implying that the basic geometric series plays a 
key role in describing relevant natural phenomena as already 
inferred by Gompertz lfl2ll . In other words, the root causes 
do not add but multiply among them in many systems, in- 
cluding economic and social ones. Our analysis accounts for 
the stochastic origin of ubiquity of Gompertzian growths, sug- 
gesting also a method to determine the order of magnitude of 
the microscopic stochastic parameters, and in particular of the 
diffusion parameter, by the statistical observation of macro- 
scopic sizes. 

Some systems, as for example cell aggregates, develop by 
a "birth-death" process, i. e. more general by the competition 
between aggregation and disaggregation. For this system, the 
basic process is then a branching process leading in a suitable 
limit to a log-normal behavior. Some of these aggregates are 
characterized by self-regulation: an example is that of tumor 
cells forced to grow in three dimensions, which cannot go be- 
yond a critical diameter regardless of how often new medium 
is provided or how much open space is made available. These 
growths cannot thus be described by a stochastic background 
where the parameters are a priori defined and one must resort 
to a dynamical setting. In the last section of this work we 
then propose a dynamical conceptual framework that include 
self regulation mechanisms. We exploit a stochastic varia- 
tional principle, whose canonical structure generates a control 
equation ruling the dynamical update of the forward velocity. 
Self regulation is then implemented by choosing the poten- 
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tial function as a (nonlinear) functional of the density of the 
system, so introducing a self consistent feedback of the rela- 
tive density variations. This conceptual framework opens the 
way to a reliable control of the growing system by external 
actions which, for example, suitably modify, on some charac- 



teristic time scales, the density profile of the system. One can 
thus effectively develop practical methods to reverse the grow- 
ing trend, reducing, for example, a cancer size. The general 
scheme here outlined will be applied and tested in forthcom- 
ing papers, focusing on the behavior of specific systems. 
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